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ABSTRACT 

We report three RXTE/PCA observations of the low-mass X-ray binary 4U 0614+09. 
They show strong (~30% rms) band-limited noise with a cut-off frequency varying 
between 0.7 and 15 Hz in correlation with the X-ray flux, f x . We observe two non- 
simultaneous 11-15% (rms) kHz peaks near 728 and 629 Hz in the power spectra of two 
of our observations when / x ~10 _9 erg cm~ 2 s _1 (2-10 keV), but find no QPO (<6% 
rms) when f x is half that. We suggest that count rate may not be a good measure for 
M even in sources as intrinsically weak as 4U 0614+09, and that QPO frequency and 
noise cutoff frequency track M more closely than count rate. The QPO increases in 
rms amplitude from (11+1.3)% at 3 to (37+12)% at 23 keV; the fractional amplitude 
of the band-limited noise is energy- independent. This suggests different sites of origin 
for these two phenomena. The spectrum of the oscillating flux roughly corresponds to a 
black body with temperature (1.56+0.2) keV and radius (500+200) m (other models fit 
as well) , which might indicate the oscillations originate at a small region on the neutron 
star surface. 
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1. Introduction 

Recently, kilohertz quasi-periodic oscillations (kHz 
QPOs) have been discover ed in e leven low-mass X-ray 
binaries (see van der Klis 1997 for a review). Often, 
the X-ray power spectra show twin kHz peaks mov- 
ing up and down in frequency together. Sometimes a 
third kHz peak is detected near a frequency equal to 
the separation frequency of the twin peaks, or twice 
that, suggesting a beat-frequency interpretation, with 
the third peak near the neutron star spin frequency 
(or twice that). However, in Sco X-l the twin peak 
separation varies, which is not consistent with a sim- 
ple beat-frequency interpretation (van der Klis et al. 



1997) 



al. 

rao 

al. 



1978; ; Brandt et 



In the X-ray burster (Swank et al. 
1992|) and suspected atoll source (Singh & Appa- 
1995[ ) 4U 0614+09 twin kHz peaks occur (Ford et 
1997 ). The peaks move between 480 and 800 Hz 
and 520 and 1150 Hz, respectively. Their separa- 
tion is consistent with being constant near 323 Hz. 
There has been one 3.3a detection of a third peak 
near 328 Hz (in the persistent emission, whereas in 
the three other sources where a third peak occured it 
was seen during X-ray bursts). 

It is yet unclear how the properties of the kHz QPO 
relate to type and state of the sources in which they 
have been observed. In this paper we analyze new X- 
ray timing and spectral data on 4U 0614+09 (prelim- 
inary reported previously in van der Klis et al. |1996 ). 
We point out a number of correlations between the 
QPO properties and those of the broad-band noise 
and the X-ray spectra, and present for the first time 
an analysis of the photon energy spectrum of the os- 
cillating flux in kHz QPO. 

2. Observations 

We observed 4U 0614+09 with the proportional 
counter array (PCA) on board NASA's Rossi X-Ray 
Timing Explorer (Bradt, Rothschild, & Swank |l993p 
three times, in February, March, and April 1996 (Ta- 
ble 0). We simultaneously collected 2-60 keV data 
with a time resolution of 8 /is in 8 energy bands, and 
16 s in 129 bands. The background and deadtime 
corrected count rates were 235, 571, and 256 c/s, re- 
spectively. A change in PCA gain between March and 
April affected these values only slightly. In April, only 
3 out of the 5 PCA detectors were active; the 5 de- 
tector count rate was 426 c/s. We observed no X-ray 
bursts. 



We calculated power spectra of the 8 /xs data and 
subtracted the Poisson noise and the Very Large 
Event window contribution (Zhang et al. 1995; Zhang 
1995| ). We obtained background measurements from 
slew and Earth occultations data, and used them 
to renormalize all power spectra to fractional rms 
squared per Hertz (see van der Klis |l995| ), and to cor- 
rect the X-ray spectra obtained from the 16 s data. 

3. Results 

All power spectra (Fig. |l|) show strong (^30% rms) 
band-limited noise, which we fitted with a broken 
power law. In two cases we observed kHz QPOs. 
These we fitted with Lorentzian peaks. The fit results 
are listed in Table Q. The power-law break frequency 
varied from 0.7 Hz in February, via 15.4 Hz in March 
to 6.6 Hz in April. 

During March and April strong QPO are present 
near 727 and 629 Hz, respectively. The QPO proper- 
ties did not vary significantly within each observation. 
We detect no other kHz QPO peaks. In particular, 
any peaks 323 Hz above or below our detected peaks 
are 3 to 14 times weaker than these (Table [l]). The 
upper limit on any 328 Hz peak is 3.5% (rms; 95% 
confidence) . 

Countrates in March were sufficient to study the 
photon energy dependence of the QPO. We fitted 
power spectra in 5 energy bands (1.1-4.3-6.8-8.4- 
15.5-69.8 keV) with the break frequency, the power 
law slopes, and QPO frequency and width fixed to 
the 2-60 keV values (none of these varied significantly 
with photon energy). As shown in Fig. 2a the rms am- 
plitude of the QPO rose from 11% at 3.3 keV via 20% 
at 11 keV to (37+12)% at 22.8 keV. Similar energy 
dependencies were seen in other sources. The frac- 
tional rms amplitude of the band-limited noise did 
not vary significantly as a function of photon energy. 

Fig. 2b shows the energy spectrum of the oscillat- 
ing flux (QPO rms amplitude in units of c/s/keV vs. 
photon energy). For reference, we quote the results 
of a blackbody fit: the best fit (x 2 = 14.4 with 3 
dof) has a temperature of (1.56+0.2 keV) an d a ra - 
dius of (500+200 m) (at 3 kpc, Brandt et al. |l992| ). 
The resolution of the spectrum is low and many spec- 
tral models are consistent with it. For example, the 
data can be fitted by a ~ 2.5% temperature variation 
in a ~ 1.1 keV blackbody spectrum with a radius of 
10 km, or (x 2 = 4.9 with 4 dof) by ~ 5% optical 
depth variations in an unsaturated Comptonization 
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spectrum. (This latter result was obtained by fitting 
the fractional rms spectrum.) 

To convert count rates to fluxes, we fitted the 2- 
50 keV energy spectra with a blackbody plus a power 
law modified by interstellar absorption (Table [l]) . In 
April and during a dip in soft color in March (below), 
a power law alone could fit the spectrum. The inferred 
2-10 keV fluxes were 0.4, 1.0, and 0.75 xl0~ 9 erg 
cm~ 2 s _1 in February, March, and April, respectively; 
the 2-50 keV fluxes are 50, 30, and 40% higher, re- 
spectively. As TVh is low absorbed and "unabsorbed" 
fluxes are the same. We tried various other spectral 
shapes, with no effect on the derived fluxes. Count 
rate, flux nor power spectra were affected by the soft 
color dip. 

To compare 4U 0614+09 to confirmed atoll sources 



(cf. Hasinger & van der Klis 198E ), we produced an X- 



ray color-color diagram using the three spectral bands 
2-5-10-50 keV (Fig |). We corrected for detector 
gain changes by comparing the count rate ratios with 
incident flux ratios obtained in the same bands from 
the spectral fits. Once again we tried different spec- 
tral shapes and found negligible effects on the color 
corrections. There was a ~ 800 s dip in soft color 
starting 5200 s into the March observation. Apart 
from this, the colors within each observation show no 
significant changes, with upper limits on any intrin- 
sic color variations of ~10%. Except for the March 
excursion, this is entirely consistent with atoll-source 
island-state behavior. 

4. Discussion 



From color-color diagrams Singh & Apparao (|l995| ) 
suggested 4U 0614+09 to be an atoll source. Our 
color-color diagrams and spectral fits, and simultane- 
ous strong band- limited noise with 0.7-15 Hz cutoff 
frequencies confirm that the correlated X-ray spec- 
tral and timing properties of 4U 0614+09 are those 
of an atoll source, in the island state during our ob- 
servations. The large fractional amplitude and low 
cut-off frequency of the band-limited noise and the 
relatively hard X-ray spectra make the source simi- 
lar to the atoll sources 4U 1608-52 and 4U 17 05-44 
in the island state (Hasinger & van der Klis |l989j ; 
Langmeier, Hasinger, & Trimmer 1989; Yoshida et 
al. 1993; Berger & van der Klis 1996| ), which in 
turn resemble the black hole candidates Cyg X-l and 
GX 339-4 in the low and intermediate states (van der 
Klis |!994a|; Berger & van der Klis |1996|; Belloni et al. 



1996| ; Crary et al. [1996| ; Mendez & van der Klis |1997| ). 
All these sources are quite hard. Their energy spectra 
fit a soft component with kT ~ 1 keV plus a power 
law with photon index 1.6-2.5. Their power spectra 
show band limited noise with a ~ 0.1 — 10 Hz cut-off 
frequency that is anticorrelated with the level of the 
flat top (Belloni & Hasinger 1990; Mendez & van der 
Klis 1997). This anti-correlation holds also for our 



power spectra of 4U 0614+09. Break frequencies and 
fractional rms at the break are fully consistent with 
the existing relation between these quantities from 
other sources (cf. Fig. 3 of Mendez & van der Klis 



1997|) 



If as proposed (van der Klis 1994b; Mendez & van 



der Klis 1997), the break frequency of the band lim- 
ited noise component here is an indication for Al , then 
M increased from February to March, and then de- 
creased to an intermediate value in April. This is 
consistent with the variations in X-ray flux among 
our three observations. Note, that this is not always 
the ca se in similar sources (e.g., 4U 1608-52; Yoshida 
et al. |1993D . 

In this picture, we observe no kHz QPO (<6% 
rms), when M is lowest, and 11-15% amplitude QPOs 
at higher inferred M. Although kHz QPOs often get 
stronger when inferred M drops (in 4U 1636-53, Wi- 
jnands e t al. |1997|; KS 1731-260, Wijnands & van 



der K lis |1997| ; 4U 1820-30, Smale, Zhang, & White 
1997), apparently in 4U 0614+09 there is a value of 



M below which the QPO becomes weaker. 

We now turn to the identification of our QPO 
peaks. In our data at most one peak is present at 
each time, whereas for similar count rates Ford et 
al. (1997) usually find twin peaks. The count rate 
vs. QPO frequency relations do not clearly identify 
our peaks as either the higher- or the lower-frequency 
ones (Fig. ^). The QPOs amplitudes do give a clue. 
Ford et al. (1997) found the lower- frequency peak to 
have an rms amplitude 0.25-0.75 times that of the 
higher-frequency peak when the count rate was near 
400 c/s, and 0.5-1.5 that at 600-700 c/s. If our peaks 
are the higher-frequency ones, these ratios are <0.3 
to < 0.4 in our data. If they are the lower-frequency 
ones, these values are >2.4 to >3.5. As our countrates 
are 400-600 c/s, our peaks are probably the higher- 
frequency ones. Then, in our data we have a positive 
correlation between QPO frequency and count rate, 
but one that does not fit either of the relations ob- 
served previously (see Fig. ||). 

The QPO frequency vs. count rate diagram of 
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4U 0614+09 (Fig. ||) shows that there is no single 
relation that describes the correlation between those 
two quantities. This means that either countrate or 
frequency (or both) does not track M well. We sug- 
gest that count rate is not a good measure of M. Both 
spectral changes and variations in the anisotropy of 
the emission can destroy the expected correspondence 
between these two quantities. To the extent that 
bolometric corrections can be accurately performed, 
a conversion to energy flux adjusts for the spectral 
changes. However, the experience from the Z sources 
(at admittedly much higher M) shows that even de- 
rived bolometric X-ray fluxes do not always track M 
well. In 4U 1608-52, in a very similar series of island 
state observatio ns as the present ones of 4U 0614+09, 
Yoshida et al. ( 1993| ) found that the bolometric flux 
did not exhibit a one-to-one correspondence to cutoff 
frequency of the band-limited noise. Whereas there 
are many processes that can quite easily change the 
observed X-ray count rates, colours, spectral param- 
eters, and even bolometric fluxes from their original 
values, the noise cutoff frequency is not so easily af- 
fected. The noise cutoff can, in principle, be lowered 
by scattering delays. The kHz QPO frequency, on 
the other hand, is very likely a direct diagnostic of 
the dynamics of the inner flow and is therefore very 
hard to modify by any propagation effect. It is pos- 
sible, then, that QPO frequency is well correlated to 
M, but countrate is not. It will be of great interest to 
see if the correlation suggested by the joint decrease of 
the noise cutoff frequency and QPO frequency from 
March to April (see Table [l]) will hold up in future 
analyses of other data, or whether a QPO frequency 
vs. flux relation will turn out to be the more repro- 
ducible, as this may provide clues to both the best 
measure of Al in atoll sources in extreme island states 
(and perhaps black-hole candidates) and to the phys- 
ical origin of kHz QPO and broad-band noise in these 
systems. 

We observed a very strong energy dependence in 
the kHz QPO but no energy dependence in the band- 
limited noise. One plausible interpretation of this, 
namely, that the low energy photons undergo scatter- 
ing with a characteristic delay time scale intermedi- 
ate between the QPO time scale (1 ms) and that of 
the noise (0.1-10 s), can probably be excluded on the 
basis of the extremely small (20 fxs and <50 /is) time 
lags between the kHz QPO signals at different photon 
energies recently reported by Vaughan et al. ( |1997 ). 
(Of course, scattering with shorter characteristic time 



scales cannot be excluded and is in fact likely on other 
grounds.) The kHz QPO spectrum which resembles 
a blackbody shape with a characteristic temperature 
of 1.6 keV and radius of 500 m, might indicate an 
origin associated with a relatively small area on the 
neutron star surface, whereas the band-limited noise 
may have a different site of origin, perhaps in the 
inner disk. This would be in accordance with the ob- 
servation that band-limited noise is a common trait 
among neutron stars and black holes (van der Klis 
1994a ), whereas correlated twin kHz QPO peaks are 
so far unique to neutron star systems. 

This work was supported in part by the Nether- 
lands Organization for Scientific Research (NWO) un- 
der grant PCS 78-277 and by the Netherlands Foun- 
dation for research in astronomy (ASTRON) under 
grant 781-76-017. MM is a fellow of the Consejo Na- 
tional de Investigaciones Cientificas y Tecnicas de la 
Rcpublica Argentina. WHGL acknowledges support 
from the National Aeronautics and Space Adminis- 
tration. JVP acknowledges support from the Na- 
tional Aeronautics and Space Administration through 
contract NAG5-3269. FKL acknowledges support 
from NSF, through grant AST 93-15133, and NASA, 
through grant 5-2925. 



4 



REFERENCES 

Belloni, T., & Hasingcr, G. 1990, A&A, 227, L33 

Bclloni, T., Mcndez, M., van der Klis, M., Hasinger, 
G., Lcwin, W.H.G., & van Paradijs, J. 1996, ApJ, 
472, L107 

Bcrger, M., & van der Klis, M. 1996, A&A, in press 

Bradt, H. V., Rothschild, R. E., & Swank, J. 1993, 
A&AS, 97, 355 

Brandt, S., Castro-Tirado, A. J., Lund, N., Dremin, 
V., Lapshov, L, & Sunyaev, R. 1992, A&A, 262, 
L15 

Crary, D. J., et al. 1996, ApJ, 462, L71 

Cui, W., Hcindl, W. A., Rothschild, R. E., Zhang, S. 
N., Jahoda, K., & Focke, W. 1997, ApJ, 474, L57 

Ford, E., et al. 1997, ApJ, 475, L123 

Hasingcr, G., & van der Klis, M. 1989, A&A, 225, 79 

Langmeier, A., Hasinger, G., & Triimper, J. 1989, 
ApJ, 340, L21 

Mendez, M., & van der Klis, M. 1997, ApJ, 479, 926 

Singh, K. P., & Apparao, K. M. V. 1994, ApJ, 431, 
826 

Smalc, A. P., Zhang, W., & White, N. E. 1997 
IAU Circ.6507 

Swank, J. H., Becker, R. H., Boldt, E. A., Holt, S. S., 
& Serlemitsos, P. J. 1978, MNRAS, 182, 349 



Vaughan, B. A., et al. 1997, ApJ, in press 

Wijnands, R. A. D., & van der Klis, M. 1997, ApJ, 
482, L65 

Wijnands, R. A. D., van der Klis, M., van Paradijs, 
J., Lcwin, W. H. G., Lamb, F. K., Vaughan, B. A., 
& Kuulkers, E. 1997 ApJ, submitted 

Yoshida, K., Mitsuda, K., Ebisawa, K., Ucda, Y., Fu- 
jimoto, R., Yaqoob, T., & Done, C. 1993, PASJ, 
45, 605 

Zhang, W. 1995, XTE/PCA Internal Memo, 5-23-95 

Zhang, W., Jahoda, K., Swank, J.H., Morgan, E. H., 
& Giles, A. B. 1995, ApJ, 449, 930 



van der Klis, M. 1994a, ApJS, 92, 511 

van der Klis, M. 1994b, A&A, 283, 469 

van der Klis, M. 1995, in X-ray Binaries, eds. W.H.G. 
Lewin, J. van Paradijs & E.P.J, van den Heuvel, 
(Cambridge U. Press, Cambridge) p. 252 

van der Klis, M., van Paradijs, J., Lewin, W. H. G., 
Lamb, F. K., Vaughan, B. A., Kuulkers, E., & Au- 
gusteijn, T. 1996, IAU Circ.6428 

van der Klis, M. 1997, in the Proceedings of the Wise 
Observatory 25th Anniversary Symposium "Astro- 
nomical Time Series" , Tel Aviv, in press 

van der Klis, M., Wijnands, R. A. D., Home, K., & This 2 - C olumn preprint was prepared with the AAS IATeX 

Chen, W. 1997, ApJ, 481, L97 macros v4.0. 



5 



Table 1 

Power- and X-ray spectral parameters 



1996 Feb 26 1996 Mar 16 1996 Mar 16 1996 Apr 13 

(soft color dip) a (main) a 



Start time (UTC) 


08:53 


09:53 


08:26 


09:50 


Exposure time (ks) 


5.0 


0.8 


4.3 


5.2 


Power spectra 


BLN rms [%] b 


30.2 ±0.9 


29.4 ± 1.6 


28.1 ± 0.7 


35.3 ± 1.1 


fbrcak [Hz] 


0.7 ±0.1 


15.4 ± 1.5 


15.4 ± 0.8 


6.6 ±0.5 


Pbrcak rms [%]° 


12.2 ± 1.1 


3.0 ± 0.2 


3.2 ± 0.2 


4.5 ±0.2 


QPO rms [%] 


< 6 


11.0 ± 1.8 


13.9 ± 0.7 


14.6 ±0.5 


Freq. [Hz] 




728 ± 8 


719 ± 5 


629 ± 4 


FWHM [Hz] 




54 ± 27 


99 ± 14 


98 ±9 


2nd. QPO rms [%] 




< 4 


< 4 


< 6 


X 2 /dof d 


284/229 


270/346 


329/346 


259/242 


Energy Spectra 


N H [10 _22 cm- 2 ]° 


(< 1.34) 


(< 0.22) 


(< 0.10) 


0.14 (< 0.39) 


kT [kcV] 


0.42 ± 0.05 




1.10 ± 0.18 




n (photon-index) 


1.97 ± 0.03 


2.54 ± 0.03 


2.41 ± 0.03 


2.30 ± 0.03 


PL Flux f 


0.36 ± 0.02 


1.00 ± 0.05 


0.89 ± 0.08 


0.75 ± 0.02 


BB Flux f 


0.04 ± 0.01 


<0.08 


0.12 ± 0.05 


<0.02 


X 2 /dof d 


103/90 


92/94 


80/92 


95/83 



a The March observation was divided into 2 parts; see text for details. 
b Band limited noise. 

c Powcr of the BLN component at inbreak- 
d Numbcr of degrees of freedom. 

°Nh was consistent with zero in all of our fits. The best-fit values arc given first, followed by 
the upper limits in parentheses. 

f 2— 10 keV flux in units of 10 -9 ergs cm -2 s — 1 . When only an upper limit is given for the 
blackbody flux, this parameter was kept fixed at zero in the final fit. 

Quoted errors represent 90% confidence intervals for the fits to the X-ray spectra and la con- 
fidence intervals for the fits to the power spectra. Quoted upper limits are 95% confidence. A 2% 
systematic uncertainty was included in the X-ray spectral errors to account for the calibration 
uncertainties (Cui et al. 0.9971). 
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Fig. 1. — The power spectra for the 3 observations 
mentioned in the paper. The data from February, 
March, and April are represented by open circles, 
filled circles, and squares respectively. The signifi- 
cance of the peak at ~ 200 Hz in the power spec- 
trum from April is less than 3a and was not fitted. 
We included a QPO at 6.5 Hz in fitting the February 
spectrum. 
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Fig. 2. — (aj Fractional rms amplitude vs. photon 
energy spectrum of the QPO from March 1996. (b) 
Energy spectrum of the oscillating flux, with black- 
body fit, for the same data set. 
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Fig. 3. — Color-color diagram showing all the ob- 
servations mentioned in the paper. Symbols are the 
same as in Fig. ^|. Each point represents 16 s of data. 
Typical error bars for each observation are indicated. 
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Fig. 4. — The count rate vs. QPO frequency diagram 
of all 4U 0614+09 data reported to date. Pluses, as- 
teriks, and po wer-l aw relations have been taken from 
Ford et al. ( 1997 ); filled circles represent the new 
data presented in this paper. 
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